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Abstract. The detection of high energy nentrinos (10^^ — 10^*^ eV) is an important step 
toward nnderstanding the most energetic cosmic accelerators and wonld enable tests of fnn- 
damental physics at energy scales that cannot easily be achieved on Earth. In this energy 
range, there are two expected popnlations of nentrinos: the astrophysical flnx observed with 
IceCnbe at lower energies (~ 1 PeV) and the predicted cosmogenic flnx at higher energies 
(~ 10^® eV). Radio detector arrays snch as RICE, ANITA, ARA, and ARIANNA exploit the 
Askaryan effect and the radio transparency of glacial ice, which together enable enormons 
volnmes of ice to be monitored with sparse instrnmentation. We describe here the design for 
a phased radio array that wonld lower the energy threshold of radio techniqnes to the PeV 
scale, allowing measnrement of the astrophysical flnx observed with IceCnbe over an extended 
energy range. Meaningfnl energy overlap with optical Cherenkov telescopes conld be nsed for 
energy calibration. The phased radio array design wonld also provide more efficient coverage 
of the large effective volnme reqnired to discover cosmogenic nentrinos. 
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1 Introduction 

The IceCube experiment has recently detected neutrinos of astrophysical origin with energies 
up to ~2 PeV [1]. This is a breakthrough in high energy astrophysics, and represents the 
first observation of astrophysical neutrinos (except those seen from the sun and supernova 
1987a [2, 3]). There are several candidate source classes that could make neutrinos up to 100 
PeV, including active galactic nuclei, gamma-ray bursts, and starburst galaxies (for reviews, 
see [4-6]). In the coming years, we seek to identify the currently unknown sources of PeV 
astrophysical neutrinos and determine the energy spectrum of those sources to learn about 
the physics that drives their central engines. 

In addition, the detection of ultra-high energy (UHE) neutrinos (> 100 PeV) would 
open a new window into the universe, allowing us to study distant energetic astrophysical 
sources that are otherwise inaccessible. UHE neutrinos are created as a byproduct of the 
so-called GZK process [7], the interaction of a cosmic ray (E > 10^®'® eV) with a cosmic 
microwave background photon [8, 9]. The current best limit on flux of UHE neutrinos comes 
from IceCube below 10^®'^ eV [10], and ANITA above 10^®'® eV [11]. Detection of these 
neutrinos, often called cosmogenic neutrinos, would shed light on the source of the highest 
energy cosmic rays, tell us about the evolution of high energy sources in our universe, and 
give us information about the composition of high energy cosmic rays. We expect neutrinos 
from the GZK process to have extremely high energies, above 100 PeV [12-17]. Detection 
of UHE neutrinos would also allow us to study weak interactions at center of mass energies 
> 100 TeV, much greater than those available at particle colliders on Earth, such as the 
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LHC [18-20]. As extremely relativistic particles, UHE neutrinos would enable sensitive tests 
of Lorentz invariance [21] and other beyond Standard Model processes. 

There are clear requirements for an ideal high energy neutrino observatory driven by 
the twin science goals of determining the origin and spectrum of the observed astrophysical 
IceCube signal, and discovering cosmogenic neutrinos to determine the origin of the highest 
energy cosmic rays while also probing particle physics at extremely high energies. The re¬ 
quirements for a high energy neutrino observatory are to 1) have the sensitivity in the PeV 
energy range to measure the observed IceCube astrophysical neutrino spectrum and extend 
the measurement to higher energies, 2) have the geometric acceptance at extremely high ener¬ 
gies (E > 10^”^ eV) to detect cosmogenic neutrinos, even in the most pessimistic of reasonable 
neutrino flux models, 3) have the pointing resolution required (sub-degree) to determine the 
origin of the observed neutrinos, and 4) have the energy resolution (factor ~ 2) required to 
measure the neutrino energy spectrum at PeV energies and above to learn about the physics 
that drives the sources of high energy cosmic rays and neutrinos. 

One promising way to detect the highest energy neutrinos is through the coherent, 
impulsive radio emission from electromagnetic showers induced by neutrino interactions in 
a dielectric — the Askaryan effect [22]. When a neutrino interacts with a dielectric such as 
ice, an electromagnetic cascade is initiated, and a net negative charge excess develops in the 
medium. This net charge excess moving faster than the speed of light in the medium yields 
Cherenkov radiation. For long wavelength, low frequency emission (frequency <1 GHz), 
the emission is coherent, and for high energy showers, the radio emission dominates. Beam 
test measurements [23-25] confirm that the emitted radio power scales as the square of the 
particle cascade energy and validate the expected angular emission pattern and frequency 
dependence from detailed numerical simulations [26]. A good medium for detection of high 
energy neutrinos is a large volume of a dense dielectric with a long radio attenuation length, 
such as glacial ice, which has an attenuation length Lq ~ 1 km [27, 28]. 

We introduce here a new type of radio detector for high energy neutrinos that will meet 
these requirements. In Section 2, we review current radio detector techniques and describe 
the phased radio array concept. The projected gains in sensitivity using a phased array are 
presented in Section 3. We conclude in Section 4. 

2 A new experimental approach: an in-ice phased radio array 
2.1 Defining a general approach 

A high energy neutrino observatory that achieves the goals outlined in Section 1 must have 
sensitivity over a broad energy range from PeV to 10® PeV scales. To reach the lowest possible 
energy threshold, the instrument should be located as close to the neutrino interactions as 
possible. The electric field strength falls as 1/r (where r is the distance from the detector to 
the shower induced by the neutrino interaction) and also suffers attenuation in the detection 
medium. The lowest threshold will be achieved by a set of detectors that is directly embedded 
in a detection medium with a long radio attenuation length, such as glacial ice. 

To achieve both a low energy threshold and increase the effective volume, the signals 
used to trigger the detector should have as high an effective gain as possible. High-gain 
antennas produce a higher signal-to-noise ratio (SNR) per antenna compared to low-gain 
antennas for radio signals aligned with the boresight, but have reduced angular coverage and 
may be impractical to deploy down a narrow borehole. Combining signals from multiple 
low-gain antennas in a phased array provides another way to achieve high gain and therefore 
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allow weaker neutrino signals to be detected while also allowing full angular coverage. The 
phased array approach is the topic of this work. 

2.2 Current radio detector techniques 

RICE was a pathfinder experiment for the radio detection of UHE neutrinos and demonstrated 
the feasibility of many operation-critical technologies [29]. The extremely high energy range, 
above 10^® eV, is currently probed by the ANITA high altitude balloon experiment [30]. The 
proposed balloon-borne EVA experiment [31] would be a novel way to reach the highest energy 
neutrinos, above 10^® eV. The ARA [32] and ARIANNA [33] experiments, ground-based radio 
arrays in early stages of development with a small number of stations deployed in Antarctica, 
have energy thresholds > 50 PeV, reaching the heart of the cosmogenic neutrino regime. This 
is achieved by placing the detectors in the ice. 

ANITA, ARA, and ARIANNA all use a similar fundamental experimental design: an 
array of antennas (16 in the case of ARA) and a data acquisition system comprise a single 
station. The stations are quasi-independent in that each individual station can reconstruct 
a neutrino event. ANITA can be viewed as a single-station experiment in this context. For 
ground-based experiments, multiple stations can be positioned several kilometers apart to 
cover large volumes of ice, and the neutrino event rate increases linearly with the number of 
stations. 

In current and previously-deployed experiments, a threshold-crossing trigger is used 
to determine when individual antennas receive an excess in power above typical thermal 
noise [30, 34, 35]. If a sufficient number of coincident antenna-level triggers occur within 
a short time window, a station-level trigger is formed, and the antenna waveforms of the 
candidate neutrino event are digitized and recorded. Essentially this type of combinatoric 
threshold-crossing trigger is only sensitive to the amount of power seen by individual antennas 
as a function of time. With this triggering approach, the smallest signal that a station can 
see is determined by the gain of each individual antenna (i.e., how much power the antenna 
sees from signal compared to thermal noise). 

2.3 An in-ice phased array concept 

We present here a new concept for radio detection of high energy neutrinos that will for the 
first time have sensitivity in the 1 PeV energy range, allowing us to characterize the measured 
IceCube astrophysical neutrino spectrum, extend the measurement to higher energies, and 
achieve meaningful overlap with IceCube in energy for energy calibration. This new radio 
detector would also achieve improved sensitivity in the UHE regime per station, provide supe¬ 
rior pointing resolution at all energies, and provide stronger rejection against anthropogenic 
radio frequency interference compared to currently-implemented radio techniques. 

The key to lowering the energy threshold of a radio experiment and increasing sensitivity 
at higher energies is the ability to distinguish weak neutrino-induced impulsive signals from 
thermal noise. For antennas triggering independently, the amplitude of the thermal noise is 
solely determined by the temperature of the ice and the noise temperature of the amplifiers 
(a smaller effect than the ice itself). 

Combining signals from many antennas in a phased array configuration averages down 
the uncorrelated thermal noise from each antenna while maintaining the same signal strength 
for real plane-wave signals (such as neutrinos). If we combine the signals from multiple 
antennas with the proper time delays to account for the distance between antennas, we 
can effectively increase the gain of the system of antennas for incoming plane waves from a 
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given direction. Many different sets of delays with the same antennas can create multiple 
effective antenna beam patterns that would together cover the same solid angle as each 
individual antenna but with much higher gain. This procedure is called beamforming, and is 
an economical and efficient way to achieve the extremely high effective gain needed to push 
the energy threshold down to 1 PeV. Such interferometric techniques have been extensively 
used in radio astronomy (for a review, see [36]). 

The effective gain Geg in dBi or dBd of the system is determined by the gain of each 
individual antenna in dBi or dBd, G, and the number of antennas, N, by: 

Geff = 10 logio(Nx 10^/10). (2.1) 

The trigger threshold of radio detectors is typically set by the rate at which antenna 
waveforms can be digitized and recorded while maintaining a high livetime fraction. Lower 
thresholds on the electric field at the antennas correspond to both increased efficiency for 
neutrino signals and higher trigger rates. For a phased array, each trigger channel corresponds 
to a single effective beam, and the station-level trigger is the union of simple threshold-crossing 
triggers on the individual effective beams of the phased array, rather than individual antennas. 
In this configuration, the threshold on the electric field could be reduced by roughly the square 
root of the number of antennas that are phased together while maintaining the same overall 
trigger rate per trigger channel. Since the electric field produced at the antenna from Askaryan 
emission scales linearly with the energy in the particle cascade, this reduction in the effective 
electric field threshold directly translates into a lower energy threshold for finding neutrinos. 
Equivalently, the effective volume of the detector is increased at fixed neutrino energy since 
events could be detected from farther away. To minimize the number of trigger channels, 
thus minimizing complexity and cost, the phased array that provides the trigger needs to be 
as closely packed as possible. Since the overall physical size of the array compared to the 
wavelength of radiation determines the angular resolution of the instrument, the closer the 
spacing between antennas, the larger each effective beam is for a given frequency of radiation, 
and the fewer channels are required to cover the same solid angle of ice [36]. 

Further sensitivity gains may be possible by recognizing that the antennas that form 
the trigger do not have to be the same antennas used for detailed event analysis. Indeed, it 
is advantageous to construct two distinct, co-located antenna arrays. The first is the phased 
array that is as closely packed as possible and provides the most sensitive trigger possible 
(the “trigger array”). The second array is a set of antennas that are spaced as far apart as is 
reasonably possible (many tens of meters) to provide the best pointing resolution and energy 
resolution possible for neutrino events (the “pointing array”). Although the station triggers on 
the compact trigger array, we do not need to digitize and record signals from those antennas 
- instead, we digitize and record signals from the much farther spaced pointing array based 
on the timing given by the trigger. The compact trigger array also has the benefit that the 
position of the antennas only needs to be determined to within a few inches in order to phase 
the antennas properly. The distinct sets of antennas that comprise the pointing array and 
the trigger array enable both optimal sensitivity and optimal directional pointing and energy 
reconstruction compared to the scenario where one compromises both technical goals in order 
to find middle ground that uses the same antennas for the trigger and pointing. This two- 
component station design represents a departure from all previous radio detector neutrino 
experiments. 

The phased trigger array also has the benefit that it can help reject anthropogenic radio 
frequency interference. A low-gain antenna has a very wide beam pattern, and is sensitive to 
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many incoming directions of radiation. However, a high-gain antenna, such as is effectively 
created by the phased trigger array, is highly directional. Since this conhguration would 
have many effective high-gain beams, each pointed in a different direction, at any given time 
we can easily mask out directions where there is man-made interference from the trigger. 
Since man-made noise tends to come from specihc incident directions corresponding to the 
location of sources of noise, the ability to mask out directions from the trigger improves 
background rejection and allows the threshold of the trigger to be set by thermal noise for 
a larger fraction of time, rather than the threshold increasing when more man-made radio 
frequency interference is present. 

Equation 2.1 holds for the case where the signal is amplihed with a low-noise ampliher 
for each antenna before the beams are formed (i.e., the noise contribution from each channel 
is the same and no signihcant noise is introduced after the beams are formed). If the loss 
in the system is low enough, the hrst stage of amplihcation could be performed after the 
beams have been formed, so each beam has a noise contribution from only one ampliher. 
This improvement would come at the cost of including any noise from loss in the cables and 
the beamforming components. In a typical case, the array views 250 K ice and the front- 
end system has a temperature of 70 K. Even if the beamforming hardware were lossless, 
the improvement in overall system temperature is small, approaching 10% improvement only 
after phasing 400 antennas, since the noise from the ice dominates the noise from the ampliher 
system. Despite only modest potential improvements, this technique may ultimately allow 
for fewer front-end amplihers to be used, which would reduce the cost of the array. 

An interferometric technique has been previously developed and used in post-processing 
analysis for fast, impulsive radio signals from neutrino interactions. This interferometric 
technique was hrst developed for and applied to the ANITA experiment [37], and has since 
been applied to other experiments [34, 38]. Interferometric techniques have also been used 
to search for extremely high energy neutrinos, above 10^^ eV, using the Westerbork Radio 
Synthesis Telescope and the moon as a target [39]. 

There have been efforts directed toward reconhguring the triggering scheme of currently- 
deployed or soon-to-be deployed experiments, such as ANITA and ARA, to instead use real¬ 
time correlation triggering after 3-bit digitization of the Nyquist-sampled waveforms [40]. 
This technique is under development, and is expected to be used in the upcoming fourth 
hight of ANITA. There are two reasons why one should consider doing the beamforming in 
hardware. The hrst is that it leaves open the possibility of doing the signal amplihcation after 
the beams are formed to reduce the noise contribution from the front-end electronics. The 
second is that information is lost in the 3-bit digitization that is preserved by doing the full 
phasing in hardware. 

As described above, the optimization of the geometry of a station should be different for 
a phased radio array compared to currently-deployed instruments to fully exploit the power of 
the beamforming technique. We explore here two scenarios: the hrst is a 16-channel station 
much like the deployed stations for the ARA experiment but with a different geometry and 
with the phased array trigger implemented, and the second is a 400-channel station that 
would achieve a 1 PeV threshold. 

2.4 Example: a 16-channel station 

We introduce here the example of a 16-channel station that uses dipole-like antennas and is 
deployed down boreholes beneath the hrn layer of a deep glacier. Possible sites for such a 
station, or set of stations, include South Pole, where the array would be ~ 200 m below the 
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Figure 1. An example station layout for a 16-antenna phased trigger array and accompanying 
pointing array. 



Figure 2. An example layout of the radio frequency chain of a 16-antenna phased trigger array and 
accompanying pointing array. For simplicity, not all channel paths are depicted. 
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Elevation 


Figure 3. The beam pattern for one trigger channel at 200 MHz for the configuration shown in 
Figure 1. 


surface and the ice is ~ 2.8 km deep [41, 42], and Summit Station in Greenland, where the 
array would be ~ 100 m below the surface and the ice is ~ 3.0 km deep [43]. 

We simulate a conhguration where the antennas are deployed down boreholes below the 
hrn. Our simulations conhrm the hnding from ARA [32, 34] that due to the changing index 
of refraction in the hrn as the glacial ice transitions to snow, incident radio emission is often 
dehected downwards away from a detector on the surface and is therefore undetected (see A.l 
for details). In contrast, receivers deployed below the hrn layer are less affected by refraction 
and the increase in effective volume compared to a surface conhguration is large (factors of 3 
to 10 depending on the hrn prohle). The benehts of borehole-deployed antennas make them 
a cost effective approach. 

One possible station layout is shown in Figure 1. A trigger array is constructed of 
16 dipole antennas strung vertically down one borehole as close together as possible. This 
conhguration would naturally be sensitive to vertically-polarized signals, although one could 
combine signals from orthogonally-polarized antennas to create an unpolarized trigger. The 
pointing array would be constructed of additional antennas, with both horizontal and verti¬ 
cal polarization sensitivity, and would require at least two additional boreholes to uniquely 
determine the incident direction, timing, and polarization of the radio emission and thus 
the incoming direction of the neutrino. We would only need to digitize the signals from the 
pointing array antennas. An example of the radio frequency chain for a 16-antenna station 
is shown in Figure 2. 

The effective gain of such a 16-channel trigger array of dipole antennas calculated using 
Equation 2.1 is 14.2 dBi (compared to 2.15 dBi for a single dipole), which corresponds to a 
factor of ~4 in electric held threshold. The effective beam at 200 MHz for one trigger channel 
is shown in Figure 3 as a function of elevation. The FWHM of the beam of a single trigger 
channel is 5.4° in elevation with complete azimuthal coverage for this conhguration. By 
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adjusting the delays among antennas in different trigger channels, we can cover the relevant 
range of solid angle for incoming emission from visible neutrino events with only ~15 beams. 
Antennas with a moderately higher gain that still cover the relevant range of solid angles 
would lead to a higher effective gain of the system. 

2.5 Achieving An Energy Threshold of 1 PeV 

The phased array design is scalable and can be configured to achieve an even lower energy 
threshold. For example, a phased array with 400 dipole antennas would have an effective gain 
of 28.2 dBi, and would push the electric field threshold down by a factor of ~20 compared 
to currently-implemented techniques. For large numbers of antennas, the single-borehole 
configuration for the trigger array is no longer optimal. To keep the size of the trigger array 
compact, trigger antennas would be deployed down multiple boreholes. Phased arrays with 
thousands of channels are possible. 

Since the spectrum and angular distribution of Askaryan emission is largely independent 
of the neutrino energy [44], the same experimental design principles and simulation methods 
are valid over the large range of neutrino energies we consider. 

3 Comparison with current experimental techniques using an independent 
Monte Carlo simulation 



Figure 4. Effective Volume vs. Energy for 10 stations installed 100 m below the surface at Summit 
Station, Greenland. The yellow line is for 16-channel stations with no phasing, the orange line is for 
similar stations but with phasing, and the red line is for 400-antenna phased array stations. For each 
radio array configuration, the volumetric acceptance is presented at the trigger level. Black curves 
indicate the volumetric acceptance for two difference analyses with IceCube optimized for different 
energy ranges [45, 46[. 
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Figure 5. Triggered Event Rates vs. Energy for a variety of neutrino models. Triggered event rates 
for three years of observation for 10 stations installed 100 m below the surface at Summit Station, 
Greenland. 16-channel stations with no phasing are shown in yellow, 16-channel stations with phasing 
are shown in orange, and stations with 400 phased antennas are shown in red. The top two panels 
show event rates based on two possible neutrino spectra based on the IceCube observed neutrino 
flux [1]. The top left panel is an E“^-^ power law, and the top right panel is an E“^-^ power law 
with an exponential cutoff at 1 PeV. The bottom two panels show event rates based on optimistic and 
pessimistic cosmogenic fluxes [17]. 


We have developed a Monte Carlo simulation package to quantify the acceptance of 
various radio detector configurations, and more specifically, investigate the advantages of a 
phased array design. The simulation formalism and assumed physics input are described in A, 
and validation studies are presented in B. Rather than focusing on specific antenna designs, 
signal processing chains, and analysis algorithms, we have kept the simulations general by 
defining signal detection thresholds in terms of the electric field strength arriving at the 
antennas. In this case, different station configurations correspond to different electric field 
thresholds, as described above. The parameters chosen to define the electric field threshold of 
our baseline station configuration are discussed in A.3. Depending on the particular system 
deployed, the overall results could shift up or down (by less than a factor of two), but the 
relative comparisons between configurations are valid. 

Figure 4 shows the volumetric acceptance of a 10 station detector with antennas 100 m 
below the surface at Summit Station as a function of energy for three different configurations: 
16-channel stations with no phasing (yellow), 16-channel stations with phasing (orange), and 
400-antenna phased array stations (red). See A for a description of how effective volume 
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Figure 6. Flux models used for the predictions shown in Figure 5. Shown are two possible neutrino 
spectra based on the IceCube observed neutrino flux [1]: an power law and an power law 

with an exponential cutoff at 1 PeV. Also shown are optimistic and pessimistic cosmogenic fluxes [17]. 


is calculated. We also show for comparison the acceptance corresponding to two IceCube 
analyses optimized for different energy ranges, calculated using the effective area given in [45, 
46] and the neutrino interaction cross section given in A. The IceCube curves in the hgure 
include analysis efficiency, whereas the results of our simulation do not, so the curves are 
not directly comparable. The standard method, shown with the yellow line in Figure 4, is 
comparable to the approach of the ARA experiment but with only 10 stations [32]. 

For this study, we have chosen to simulate the experiment at Summit Station in Green¬ 
land. Moving the detectors to South Pole at a depth of 200 m below the surface would 
increase the acceptance by < 20% due to a longer radio attenuation length in ice at the 
South Pole [27]. We have chosen to simulate 10 stations for this study, but we note that the 
acceptance changes linearly with the number of stations, since stations trigger independently. 

Figure 5 shows the number of events detected as a function of energy for the same 
three detector conhgurations shown in Figure 4 for a variety of astrophysical and cosmogenic 
models. We show event rates for an power law based on IceCube observations [1], an 

power law with a 1 PeV exponential cutoff [1], and optimistic and pessimistic cosmogenic 
models [17]. The flux models used for Figure 5 are shown in Figure 6. Table 1 summarizes the 
expectation values for the total number of events detected with each detector conhguration 
for each model. A harder spectrum for PeV-scale neutrinos would yield a higher event rate. 

For the most pessimistic cosmogenic neutrino flux models (not shown in Figures 5, 6, 
or Table 1) with no source evolution and a pure iron composition for UHE cosmic rays [17], 
the expected event rate with 10 stations of 400 phased antennas each is ~ 0.1 event per 
year. However, recent measurements with Pierre Auger Observatory and the Telescope Array 
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Station Configuration 

Power Law 

Power Law 
with Cutoff 

Optimistic 

Cosmogenic 

Pessimistic 

Cosmogenic 

16-antenna 

0.9 

0.0 

7.7 

2.3 

16-antenna, phased 

3.8 

0.1 

19.6 

6.0 

400-antenna, phased 

18.4 

2.2 

52.9 

15.6 


Table 1. Expectation values for the total number of triggered events in 3 years for 10 stations in 
different configurations for spectra based on IceCube observations [1] and for cosmogenic models [17]. 


disfavor a significant iron fraction in the cosmic ray composition at energies up to and even 
exceeding eV [47, 48], so we do not discuss these cosmogenic models further in this 

work. 

By phasing the antennas in a 16-antenna array, we have improved the UHE neutrino 
event rate by more than a factor of two over the non-phased case, and extended the sensi¬ 
tivity to lower energies. The 16-antenna phased configuration achieves a low enough energy 
threshold to distinguish an power law extrapolation of the observed IceCube spectrum 

from one that has a cutoff at the PeV scale. Phasing more antennas lowers the threshold even 
further, and makes marked improvements in event rates, especially at low energies (see the 
400-antenna configuration in Table 1 and Figure 5). Increasing the number of stations results 
in a linear increase in the expected number of events. 

As is evident in Figure 4, 10 stations of 16-antenna phased arrays have a larger acceptance 
than IceCube above 30 PeV, and the acceptance grows faster with energy than IceCube, so 
that by 10^® eV the acceptance is an order of magnitude more than IceCube. 10 stations of 
the 400-antenna phased arrays have a larger acceptance than IceCube above 1 PeV. 


4 Conclusions 

We have described a new concept for an in-ice phased radio array that is designed to achieve 
sensitivity to the astrophysical neutrino flux at 1 PeV and above, provide energy overlap with 
IceCube for calibration, and discover cosmogenic neutrinos in an efficient way. It is worth 
noting the scalability of the radio technique for increasing acceptance at all energies. The 
acceptance increases linearly with the number of stations, and further gains can be realized 
by phasing more antennas in each station, particularly at PeV neutrino energies. An array 
of 100 stations each with 400-antenna phased arrays could detect hundreds of neutrinos at 
PeV energies and above each year. Such an experiment would revolutionize our view of the 
high-energy universe. 
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A Simulation formalism 


We used numerical Monte Carlo simulations to evaluate the acceptance of different array 
configurations considered in this work. The individual stations of the array were assumed 
to be widely spaced such that a given neutrino interaction is unlikely to be detected with 
multiple stations. In this limit, the stations act independently and it is possible to simulate 
the response of a single station and multiply the derived acceptance by the total number of 
stations in the array to obtain the total acceptance (see, e.g., [32]). We also assumed that the 
ice properties depend only on depth such that a cylindrical coordinate system with azimuthal 
symmetry can be defined. To provide a concrete example, we used Summit Station, Greenland 
as the site of the envisioned radio array and considered stations at a depth of 100 m, i.e., just 
below the firn layer [43]. 


A.l Ray-tracing 

The simulations were performed in two stages for computational efficiency. First, we generated 
a library of ray-tracing solutions for radio waves propagating through the ice sheet. Rays 
originating at the particle cascade must intersect the station in order to be detected. This 
boundary condition implies that a complete set of viable ray-tracing solutions can be found 
by following rays launched at varying zenith angles from the station in the time-reverse sense. 
The trajectories are curved due to the varying index of refraction in different layers of the 
ice, and reflections can occur both at the ice-air interface at the top of the glacier and the 
ice-rock interface at the base of the glacier. Measurements of the firn density profile (upper 
100 m) at Summit Station show that the index of refraction at ~ 300 MHz is correlated to 
the vertical density profile determined from ice cores through the relation n = 1 + kp with 
k = 8.45 X 10“^ kg“^ m^ [49]. We applied the same formula to density profiles from GISP2 
ice core data [43] to estimate the index of refraction throughout the ice sheet. Figure 7 shows 
ray-tracing solutions for a station located at a depth of 100 m. The density profile changes 
rapidly in the firn layer of the ice sheet and therefore ray bending is expected to be most 
pronounced in the upper ~ 100 m of the ice sheet. This effect results in a “firn shadow” region 
that is inaccessible except through reflections at the bottom of the ice sheet [50]. 

At short propagation distances from the station, there are three distinct sets of con¬ 
tinuous solutions that define regions in which different classes of ray-tracing solutions exist. 
The first class involves reflections off the ice-air interface. The latter two classes involve rays 
that do not reach the surface, but which can cross over each other due to different amounts 
refraction in the firn. It is useful to distinguish between these two classes to recognize regions 
where two viable ray-tracing solutions corresponding to different zenith angles exist. The 
three solution classes are doubled when also considering rays that reflect off the bottom of 
the ice sheet, resulting in a total of six solution classes. 

At each step along a given ray-tracing solution, which are spaced from 1 ns to 5 ns in 
time depending on the differential index of refraction, we record the location, propagation 
distance, direction of propagation, and the cumulative attenuation to the station separately 
in horizontal and vertical polarization. We used the radio attenuation length measurements 
at Summit Station from [27]. Fresnel conditions are applied to compute the reflected power 
at the ice-air and ice-rock interfaces. 
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Figure 7. A ray-tracing diagram for antennas 100 m below the surface of the ice at Summit Station 
in Greenland. The left shows the full depth of the glacier (3000 m), and right is a detailed view close 
to the station. For a given zenith angle of the ray at the antenna, there are three solutions for the true 
incident radiation direction as discussed in the text. Including waves that reflect off of the ice-rock 
interface at the bottom of the glacier doubles the number of solutions. The second set of solutions 
are not shown here. The color indicates the zenith angle at which the antenna sees the ray. 


A.2 Sampling neutrino interactions and electric field production 

In the second stage of the simnlations, we randomly generated nentrino interaction geometries 
within a cylindrical target volnme centered on the station and nsed the ray-tracing library 
described above to compnte the electric field that wonld be prodnced at the antennas. Nen¬ 
trino interaction vertices were selected nniformly thronghont the target volnme and arrival 
directions were nniformly distribnted over the celestial sphere (dvr sr). For compntational 
efficiency, the size of the target volnme, 14 im) was enlarged for the highest nentrino energies 
and contracted for lower nentrino energies. In all cases, the depth was matched to the thick¬ 
ness of the ice sheet (3000 m) and the radins was chosen to be large enongh to encompass all 
detectable events. 

We focnsed on nentrino-indnced hadronic cascades when compnting the Askaryan emis¬ 
sion. Mnons and tan leptons prodnced in charged-cnrrent interactions are expected to deposit 
their energy in elongated tracks that wonld be difficnlt to detect [51], and the electromagnetic 
cascade from charged-cnrrent interactions of UHE electron nentrinos can be elongated dne to 
the Landan-Pomeranchnk-Migdal (LPM) effect [26, 52, 53], resnlting in a narrow Cherenkov 
emission cone [54]. Hadronic cascades are less affected by elongation dne to the LPM effect, 
and the width of the Cherenkov cone is therefore less dependent on nentrino energy [44, 55]. 
In addition, simnlations of the LPM effect at higher energies have been shown to prodnce 
cascades with interfering pnlse trains rather than elongated pnlses [56]. With the high gain 
provided by an in-ice phased array, these signals conld be more readily differentiated from the 
thermal noise backgronnd. Estimating the increase in acceptance dne to these signals will be 
the snbject of a fntnre stndy. Either way, the approach taken here is conservative since the 
LPM effect will only become dominant for electromagnetic showers with energies > 10^® eV. 

For each candidate nentrino event, we determined viable ray-tracing solntions by in¬ 
terpolating between rays within each respective class of solntions in the ray-tracing library 
(see above). Next, for each viable ray-tracing solntion, we nsed the analytic parametrization 
of [57] to compnte the Askaryan emission at the nniqne angle between the nentrino propaga- 
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tion direction (parallel to the cascade development) and the ray connecting the cascade to the 
station. We nsed a characteristic length scale for the shower of 1.5 m, as snggested by [57]. 
When compnting the fraction of the incident nentrino energy that is visible in the cascade, we 
nsed the analytic parametrizations for the inelasticity (Bjorken-y parameter) from [19], which 
corresponds to the energy in the hadronic cascade. The inelasticity distribntions for nentral 
cnrrent and charged cnrrent interactions of both nentrinos and anti-nentrinos are expected 
to be nearly indistingnishable in the UHE range. For simplicity, we drew valnes from the 
inelasticity distribntion for charged-cnrrent nentrino interactions. 

In the UHE regime, the nentrino-nncleon interaction cross section is snfhciently large 
that a significant fraction of npward-going nentrinos interact while passing throngh the solid 
body of the Earth and are therefore not detectable. We nsed the analytic parametrizations for 
the interaction cross section from [19] and the Preliminary Reference Earth Model (PREM) 
profile for the Earth’s density [58] to compnte the probability PEarth that nentrinos snrvive 
passage throngh the Earth as a fnnction of their energy and incoming zenith angle. The 
Earth passage probability takes into acconnt the distance traversed throngh the ice sheet 
on the way to interaction vertex. For the pnrpose of compnting the nentrino interaction 
probability, we nsed the total interaction cross section for nentrinos (as opposed to anti- 
nentrinos). In the energy range considered here, the total cross section for nentrinos is a 
factor of 1.06 to 1.25 larger than for anti-nentrinos. The cnrrent simnlations do not inclnde 
the cross section enhancement from the Glashow resonance, Vf.e —)• IT“, which wonld lead to 
a farther enhancement in acceptance for electron anti-nentrinos of energy ~ 6.3 PeV [59-61]. 

The process described above effectively samples the instrnmental sensitivity over the 
detector volnme and all possible arrival nentrino directions. A minimnm of 10® candidate 
nentrino events were processed in each of 17 logarithmic steps ranging from 10^® eV to 10^^ eV 
separated by 0.5 decades in nentrino energy. The all-sky water-eqnivalent volnmetric accep¬ 
tance (km^ sr“^) was then calcnlated for each energy as 

T^r> 47rI4im ( Pi \ t k ^\ 

^ , 1 PEarth,i ^ Pdetect,i ^ ) ■> (^-1) 

N ^ \ Pwater/ 

with i iterating over nentrino events. Pdetect.i denotes probability that the electric field at 
station resnlts in a trigger (see A. 3). pi is the ice density at the interaction vertex and pwater 
is the density of water. To compnte the event rate for varions nentrino flnx models, the 
sensitivity can be eqnivalently expressed in terms the areal acceptance (m^ sr“^) as 

AH = UH//, (A.2) 

where I is the nentrino interaction length throngh water nsing the total interaction cross 
section [62]. 

A.3 Electric field threshold for a fiducial station configuration 

The final step of evalnating the nentrino acceptance is to determine which events pass the trig¬ 
ger. For generality, we nsed a simple threshold on the electric field prodnced at the antennas 
to define a trigger criterion and we compare the acceptance of varions station confignrations 
at the trigger level. For each case that we simnlated, we assnme the same fixed trigger rate 
per channel. 

We assnme a bandwidth of 100 — 800 MHz and compnte the signal electric field strength 
by integrating the Askaryan emission over this freqnency range. As explained in detail below. 
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an electric field trigger threshold of 0.15 mV/m per antenna is eqnivalent to a 2.3a trigger 
threshold for ideal dipole antennas with the assnmed bandwidth. 

The root mean sqnare voltage received in a system (Vrms) foi' thermal noise is 


Vrms = \/kBTBZ, (A.3) 

where T is the total system temperatnre, B is the bandwidth of the system, Z is the impedance 
of the system, and /cr is Boltzmann’s constant. Assnming a system temperatnre of 70 K in 
addition to the temperatnre of the ice (250 K), a bandwidth of 700 MHz, and a 50H system 
impedance, we calcnlate that Vrms = 12.6/rV. This thermal noise contribntion sets the 
threshold of the system. 

The electric held threshold is given by 


E = Y 




^0 


(A.4) 


where V is the voltage received, is the effective area of the antenna, and Zq is the 
impedance of free space [63, 64]. is dehned as 


- 47r ’ 

where G is the antenna gain, and A is the wavelength, 
threshold is therefore given by 


(A.5) 


For a 50H system, the electric held 


E = Y 


9.73 

xVg' 


(A.6) 


Assnming a perfect dipole antenna with a gain of 2.15 dBi and a fast impnlse response, 
a reference freqnency of 200 MHz, and a trigger threshold of 2.3a above thermal noise, we 
calcnlate an electric held threshold per trigger channel of 0.15 mV/m, which is the same as 
the per-antenna threshold for the 16-antenna station conhgnration withont phasing. 

The hnal acceptance for a detector mnst also inclnde the efficiency of snccessfnlly recon- 
strncting events and discriminating signal from thermal noise and anthropogenic backgronnds. 
Previons radio detection experiments have achieved high analysis efficiencies relative to the 
trigger-level efficiency, e.g., a 96% reconstrnction efficiency for ANITA [11]. Despite a lower 
signal-to-noise per antenna expected in a phased array conhgnration, the opportnnity to 
position in-sitn stations in stable, low-noise environments as well as directional rejection of 
backgronnds (Section 2.3) conld resnlt in comparably high analysis efficiencies. 


B Simulation diagnostics 

As a basic validation of the simnlation ontpnt and to develop an intnition for the instrnment 
response, we compare the distribntions of several triggered event properties for different sta¬ 
tion conhgnrations. We investigate here two different energies (1 PeV and 1000 PeV) and the 
three different experimental conhgnrations we have discnssed in this paper (16-channels not 
phased, a 16-antenna phased array, and a 400-antenna phased array). 

Fignre 8 shows the distribntion of the distance from the antennas to the interaction 
vertex as measnred along the radio signal path for each event that satishes the trigger re- 
qnirements. As expected, lower energy events occnr closer to the detector and station conhg¬ 
nrations with lower thresholds are capable of triggering on more distant events. 
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The angle at which the detector views events relative to their neutrino propagation 
direction is shown in Figure 9. The Cherenkov angle is shown for reference as a dashed line. 
As expected, higher energy events are seen farther from the Cherenkov angle. Configurations 
with lower threshold (e.g., phased arrays) also allow for events to be seen farther from the 
Cherenkov angle where the emission is weaker. 




Distance to Interaction Vertex (in) Distance to Interaction Vertex (in) 


Figure 8. The distance from the antennas to the interaction vertex for passing events at two different 
energies (1 PeV and 1000 PeV) and the three different experimental configurations we have discussed 
(16-channels not phased, a 16-antenna phased array, and a 400-antenna phased array). Lower energy 
events occur closer to the detector, and a phased system allows more distant events to be observed. 
Note that the axis ranges are different between the two panels. 



Figure 9. The angle of observation of the radio signal compared to the incident neutrino propagation 
direction for passing events at two different energies (1 PeV and 1000 PeV) and the three different 
experimental configurations we have discussed (16-channels not phased, a 16-antenna phased array, 
and a 400-antenna phased array). The Cherenkov angle is shown with a dashed line for reference. 
Higher energy events are seen farther from the Cherenkov angle, and a phased system allows events 
to be observed farther from the Cherenkov angle. 

The distribution of incident neutrino zenith angles is shown in Figure 10. The events 
are partitioned into bins of equal solid angle. At 1000 PeV, the acceptance is dominated by 
events with arrival directions above the horizon due to effective screening by the Earth for 
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up-going neutrinos. As expected, the zenith angle distribution is broader for 1 PeV neutrinos, 
for which the interaction length is comparable to the size of the Earth. 



Figure 10. The zenith angle of the incident neutrino direction for passing events at two different 
energies (1 PeV and 1000 PeV) and the three different experimental configurations we have discussed 
(16-channels not phased, a 16-antenna phased array, and a 400-antenna phased array). The arrays 
are most sensitive to neutrinos that are slightly down-going. 
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